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Abstract 


The effects of phase fractionation on the composition of oils, condensates and gases have been studied in the 
laboratory by a differential evaporation experiment. The experiment was designed to simulate the behaviour of a sin- 
gle-charge petroleum during (upward) migration or in a reservoir undergoing uplift, scenarios thought to be typical for 
the Norwegian continental shelf. After each pressure reduction step, samples from all phases (oil, condensate, free and 
associated gas) were collected and analysed. The results show that, in general, small molecules are preferentially enri- 
ched in the gas phase, both at elevated pressure and temperature and under standard conditions (15°C, 1 bar). The 
solvent power of the gas phase for large molecules is drastically reduced at lower pressures (here below approx. 300 
bar). Carbon isotope ratios of individual compound are not affected, while variations 1n isotope ratios of fractions can 
be explained by changes in the fraction compositions. Individual geochemical parameters behave as can be expected 
based on the molecular size and structure of the compounds involved. However, a generally valid, quantitative pre- 
diction of how geochemical parameters vary as a result of phase fractionation is very difficult, as this will depend not 
only on the pressure and temperature of fractionation, but also on the composition and amounts of the oil and gas that 
are involved. The observations from this study confirm the importance of phase fractionation as an important process 
to consider when interpreting geochemical data for both condensates and oils. © 2000 Elsevier Science Ltd. All rights 
reserved. 
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1. Introduction 


Phase fractionation 1s one of the processes by which 
gas-condensates! can be generated in the subsurface, 
either during migration or in a reservoir (Zhuze and 
Ushakova, 1981; England and Mackenzie, 1989; Larter, 
1990). Petroleum is generally thought to be expelled 
from its source rock as a single-phase fluid and will 
normally migrate towards shallower depths with lower 
temperatures and pressures. At depths where the pres- 
sure drops below the bubble point of the petroleum, 
separate oil and gas phases can be formed, which can 
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' Condensate (sometimes referred to as gas-condensate) is 
that part of petroleum, which is present in the gas phase in a 
reservoir, but precipitates as a liquid during protection due to 
the reduction in temperature and pressure. 


potentially follow different migration paths and end up 
in different reservoirs. Phase fractionation can also 
occur in reservoirs due to pressure reduction (uplift or 
leakage) or the addition of (dry) gas to an oil-filled 
reservoir, as the new mixture will have a higher bubble 
point than the original oil. Subsequent remigration of 
the oil and/or gas phase can result in the presence of 
several reservoirs filled with petroleum from a common 
source but with different compositions. It is important 
to recognise the effects of these physical phase fractio- 
nation processes and to distinguish them from other 
processes controlling petroleum composition, such as 
variations in source rock type and thermal maturity 
level. Only in this way can reliable correlations between 
oils, condensates and source rocks be carried out and 
a good understanding of the petroleum geology be 
obtained. 

The changes in petroleum composition introduced 
by liquid-to-gas transitions can be studied in the 


0146-6380/00/$ - see front matter © 2000 Elsevier Science Ltd. All rights reserved. 


PII: S0146-6380(00)00128-5 


1420 G.W. van Graas et al. | Organic Geochemistry 31 (2000) 1419-1439 


laboratory under controlled temperature and pressure 
conditions. Such experiments are routinely performed in 
PVT laboratories to study the behaviour of petroleum 
reservoirs during production. These so-called duiffer- 
ential evaporation experiments (Fig. 1) start with a 
sample of single-phase reservoir fluid at the pressure and 
temperature of the reservoir. Subsequently, the pressure 
is reduced to below the bubble point giving a gas and a 
liquid phase. The gas phase is removed and transferred 
to standard temperature and pressure (15°C, 1 atm) via 
a “single flash” yielding ‘free’ gas* and condensate. 
The remaining liquid phase (at high p, T) is subjected to 
a next stage of pressure reduction, etc. If a sample of the 
liquid phase is collected, this will yield oil and associated 
gas after single flash to standard conditions. 

Several laboratory experiments with gas—oil transi- 
tions have been described in the literature (Price et al., 
1983; Thompson 1987, 1988; England et al., 1991; Lar- 
ter and Mills, 1991; Carpentier et al., 1996). Much of 
this work, particularly by Thompson, focussed on the 
light hydrocarbon fraction (C¢—Cg). Data on the Cy5+ 
fraction has been reported (e.g. Larter and Mills, 1991), 
but not very comprehensively, while the gas fraction 
(C\-Cs) has been largely ignored. The aim of our 
experiment was to complement the published work by 
studying the full range of changes occurring during 
phase fractionation. This was obtained by sampling all 
gas and liquid phases during each stage of a differential 
evaporation experiment and subjecting these samples to 
a broad range of geochemical analyses. It was our 
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intention to use the results from these experiments to set 
up general guidelines for how geochemical parameters 
may change due to phase fractionation. In contrast to 
the earlier reported laboratory studies, which were star- 
ted with “artificial” reservoir fluids, composed of an oil 
and a non-related gas, the starting liquid in our experi- 
ment consists of a naturally occurring reservoir fluid. 
The experiments were intended to simulate a typical 
scenario on the Norwegian continental shelf, 1.e. upward 
migration or reservoir uplift of a single-charge petro- 
leum, rather than e.g. gas flushing as may be more 
common elsewhere. The experiment was carried out 
under realistic pressures (corresponding to depths of 
approx. 5000 m and less) while the (isothermal) tem- 
perature conditions only represent the deeper reservoir 
interval and not the cooler, shallower layers. 


2. Experimental 


A modified differential evaporation experiment was 
set up in Statoil’s Fluid Laboratories. The starting 
volume of reservoir fluid was increased compared to the 
routine set-up in order to obtain a sample of every phase 
(oil, associated gas, condensate, free gas) from each 
stage. In routine experiments the liquid phase is not 
sampled during each stage in order to obtain a good 
mass balance. However, in our experiment small ali- 
quots of the liquid phase were collected during each 
stage in the modified experiment. Finally, the number of 
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Fig. 1. Schematic overview of differential evaporation experiment. pg = bubble point. 


* This gas phase is termed ‘“‘free gas” throughout this report in order to distinguish it from the (oil-)associated gas that is also 


obtained. 
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stages and the associated pressures were selected to give 
reasonable size samples for all phases. 

The experiment was carried out with a bottom hole 
sample (BHS), collected at reservoir conditions (494.8 
bar, 164.8°C) from a well offshore mid Norway. Under 
reservoir conditions this is a single-phase fluid with a 
gas—oil ratio (GOR) of 390 Sm? gas/Sm? oil and a bub- 
ble point pressure of 412 bar. During the experiment the 
pressure was reduced from reservoir pressure to | bar in 
six stages: 390, 340, 280, 200, 100 and 1 bar; all at 
164.8°C. In a geological setting, these pressures would 
be compatible with depths ranging from approximately 
3900 m to surface (assuming hydrostatic conditions). 
Obviously, the constant high temperature deviates from 
a normal geological situation, where temperature and 
pressure will decrease simultaneously. However, careful 
evaluation of the results strongly suggests that the con- 
stant temperature does not significantly alter the results, 
with the exception of the last stage with a pressure 
reduction to | bar (see discussion below). 

At each stage the gas phase was removed and a small 
sample of the liquid phase (approx. 3 ml, 7—-10% of the 
remaining liquid phase) was taken before the next stage 
of pressure reduction was carried out. These liquid and 
gas samples were collected at the elevated experimental 
temperature and pressure. After collection, they were 
transferred to standard conditions (1 bar, 15°C) via a 
single flash, thus yielding stabilised oil, associated gas, 
condensate and free gas. All four sample types were 
obtained from every stage, apart from the last stage (no. 
6) as a “‘single flash”’ of oil from | bar/164.8°C to | bar/ 
15°C obviously does not produce an associated gas. 
Some of the samples obtained (associated gas and con- 
densate from stage 5) were very small and the analytical 
program could only partly be carried out. 

Gas samples were analysed by gas chromatography 
for the distribution of C;-Cj 9+ components, while iso- 
topic measurements were performed on CO, and the 
C,-C; hydrocarbons. Oils and condensates were sub- 
jected to the following analytical program: whole-fluid 
GC including quantification of the light hydrocarbon 
fraction, gross composition analysis by thin layer chro- 
matography (latroscan), GC and GCMS analysis of the 
saturated and aromatic hydrocarbon fractions, and car- 
bon isotope measurements of the whole oils/condensates 
and all fractions. All analyses were performed in accor- 
dance with ““The Norwegian Industry Guide to Organic 
Geochemical Analyses” (NIGOGA, 1993). 


3. Results and discussion 


3.1. Overall mass balance and bulk fluid composition 


Based on the weights of all samples collected during 
the various stages, a mass balance has been set up for 


the experiment (Table 1). Gas—oil ratios (GOR) have 
been calculated for the gas phase (= free gas/con- 
densate), the liquid phase (= associated gas/oil), and for 
the whole system (free+associated gas/oil+conden- 
sate). 

Between 8 and 14% of the liquid was evaporated 
during each stage (Fig. 2). The bulk of this material is 
represented by the free gas (4-9%), with condensate 
(0.3—2.4%; 5.6% in Stage 6) accounting for the rest. The 
residual oil after Stage 6 represents 59% of the starting 
material (after correction for liquid samples removed 
during the experiment). During all stages a total of 
22265 ml gas (free and associated) were liberated from 
the oil. This is equivalent to an overall GOR of 404 
Sm?/Sm°, which is in good agreement with a GOR of 
390 Sm?/Sm? that was measured in a single flash 
experiment. 

GORs have been calculated for the entire system, the 
gas phase and the liquid phase (Fig. 3). GORs for the 
liquid phase gradually decrease from 376 at Stage | to 
57 at Stage 5 (no associated gas is released during Stage 
6). GORs for the gas phase are much higher and 
increase from 1840 at Stage | to 18293 at Stage 5, with a 
low value of 829 at Stage 6. The latter is due to combi- 
nation of high temperature (165°C) and low pressure (1 
bar) which causes a transfer of low-boiling (Cg—Co) 
components to the gas phase. The GOR of the entire 
system parallels the GOR of the liquid phase, reflecting 
the quantitative dominance of the latter over the gas 
phase. 

The quantitative dominance of the oil phase 
throughout the experiment, as described in the mass 
balance, makes that any changes in the oil phase during 
a single stage are expected to be small, especially with 
regard to bulk parameters (effects on the gasoline-range 
components can be significant), even though the 
released condensate may have a significantly different 
composition. However, small differences during each 
stage may amount to a significant cumulative effect on 
oil composition over the entire differential evaporation 
experiment. The composition of the associated gas can 
change significantly throughout the experiment as the 
whole fraction eventually is converted to free gas. 

A similar experiment was performed by Larter and 
Mills (1991). Their starting fluid was an artificial gas- 
condensate with a GOR of 3700 Sm?/Sm?. Although 
they do not report a mass balance, this higher initial 
GOR may give rise to quantitatively different effects 
compared to our experiment, as the relative sizes of the 
phases will be very important in determining the extent 
of any gas-liquid phase fractionation effects. 

Condensate—-gas ratios (CGR; mass of condensate/ 
mass of free gas) have been calculated for the gas phase 
of each stage (Fig. 3). They show an exponential 
decrease with decreasing pressures throughout the 
experiment before ending with a high value in Stage 6 
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Stage Pressure Material removed ( at standard pressure and temperature, | bar, 15°C) Liquid phase 
(bar) transferred to 
next stage? (g) 

Free gas Condensate Oil Associated gas 

Volume Density Weight Density Weight Density Volume Density 

(ml) (g/em3)_(g) (glom3) (g) (g/cm) (ml) (mgm?) 
Start 495 70.40 
Stage 1 390 2958 0.812 1.302 0.810 2.993 0.856 1314 0.855 61.78 
Stage 2 340 5414 0.810 1.455 0.805 2.746 0.856 871 0.872 S127 
Stage 3. 280 3421 0.819 0.507 0.801 2.552 0.856 594 0.888 44.13 
Stage 4 200 2738 0.807 0.232 0.790 3.534 0.856 412 0.931 38.00 
Stage 5 100 2352 0.838 0.099 0.770 2.720 0.856 234 1.009 31.66 
Stage6 = 1 1957 1.085 L773 O75). 27286 0.872 0 — _ 

Gas to oil ratio (GOR) (Sm?/Sm?) 
Stage Pressure Oil Associated Condensate Free Associated Free gas Association 
(bar) gas gas gas + free gas 

Vol. wt. vol. wt. vol. wt. vol. wt. Oil Condensate Oil 

m™ @ @ @ @ @® @m + condensate 
Start 495 70.40 390.4° 
Stage 1 390 52.92 45.30 19889 20.85 1.61 1.30 2958 2.94 376 1840 419 
Stage 2 340 47.94 41.04 13016 13.91 1.81 1.46 5414 5.38 272 2995 370 
Stage 3. 280 44.11 37.76 8789 9.57 0.63 0.51 3421 3.44 199 5405 273 
Stage 4 200 41.03 35.12 5320 6.07 0.29 0.23 2738 2.71 130 9323 195 
Stage 5 100 38.32 32.80 2172 2.69 0.13 0.10 2352 2.42 57 18293 118 
Stage6 = 1 31,29) 27,29 0 0 2.36 Lay 1957 2.60 0 829 58 


* Actual weights and volumes of fractions obtained during the experiment are listed in the upper part of the table. The lower part contains a mass 
balance for each stage, assuming full conversion of the liquid phase (oil and associated gas) to standard condtions. 


> Calculated from the sum of the fractions in the following stage 


© Measured in a single flash of reservoir fluid (495 bar, 165°C) to standard conditions (1 bar, 15°C). 


(temperature effect, see above). This exponential trend is 
different from the linear relationship between CGR and 
pressure that was reported by England et al. (1991). The 
difference can partly be explained by the fact that Eng- 
land and co-workers measured how much oil of a given 
composition could be dissolved by a gas of fixed com- 
position, while the composition of both liquid and gas 
phases vary continuously in our experiment. 

The whole fluid composition of all samples, expressed 
as relative concentrations of C,;—C,9+ pseudo-compo- 
nents (Table 2), is shown in Fig. 4. Oils comprise 85-— 
90% Cj9+ components and change very little during 
most of the experiment (up to Stage 5). Condensates, on 
the other hand, show a gradual change towards lower 
proportions of Cj 9+ at lower pressures, starting at 75% 
in Stage 1 and decreasing to 60% in Stage 4. This 
change reflects the lower “solvent power” of the gas 


> The composition of the condensate from the Stage 5 is 
considered to be unreliable due to the very small sample size. 


phase for large molecules at lower pressures. The 
amount of Cj;)+ components is strongly reduced to 
20% in the Stage 6 condensate. At the low pressure 
during this stage (1 bar) the geologically unrealistically 
high temperature (165°C) is above the boiling point of 
components up to Cy and will cause a massive transfer 
of the C,;—Coy fraction to the gas phase and most of the 
C;5—Co part of this material will end up in the condensate 
after cooling to 15°C (standard p, T). In comparison, at 
100 bar (Stage 5) the boiling point of propane (C3) is 
approx. 160°C, which means that compounds larger 
than propane will predominantly be present in the liquid 
phase. During Stages 1 through 5 the release of con- 
densate with a different composition does not sig- 
nificantly alter the composition of the remaining oil as 
the amount of condensate is small and it contains only a 
modest amount of C6—-Cy components. During Stage 6 a 
fairly large amount of condensate 1s generated, which, in 
combination with its large proportion of Cs6—Cy compo- 
nents, significantly affects the composition of the 
remaining oil. However, as pointed out earlier, Stage 6 
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Fig. 2. Mass balance for differential evaporation experiment. 
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Fig. 3. Gas—oil ratio (GOR) for overall system, gas and liquid phases; condensate—gas ratio (CGR) for gas phase. 


is not a geologically realistic situation (low pressure, conditions. Only during Stage 6, when the condensate 
high temperature). sample also changes significantly, the free gas is enri- 

The composition of the free gases changes relatively ched in C3-C7 components. Throughout the experiment 
little during Stages | to 5, as this will largely be deter- the associated gases become gradually enriched in C,+ 


mined by the gas-condensate equilibrium at standard components. This is due to the preferential evaporation 
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Differential evaporation experiment: bulk fluid composition, gas wetness* and Thompson indices” 


a. Associated and free gases 


Associated gas 


Free gas 


S. Flash Stage 1 Stage 2 Stage 3 Stage 4 Stage 5 Stage 1 Stage 2 Stage 3 Stage 4 Stage 5 Stage 6 


Composition (wt.% ) 


N,+CO, 125 1s 11.6 11.6 112 10.4 12.1 1222 12.0 12.6 127 9.2 
C, 44.7 44.5 42.6 40.5 36.2 28.6 49.7 50.2 49.1 50.3 46.1 2561 
C> 13.4 13.9 14.3 15.1 16.0 18.1 12.9 12,7 127 13.1 14.5 16.0 
C3 11.0 11.6 123 13.4 15.3 19.3 9.8 9.3 9.6 9.4 10.8 18.0 
Cy yee 8.1 8.8 9.5 10.9 13.5 6.5 6.1 6.4 5.8 6.7 14.6 
C;° 4.3 4.4 4.7 4.8 5.3 a 3.5 3.4 3.6 3.1 3.4 8.2 
Ce 2.4 2.3 2.3 2.2 2.3 2.0 2A 2,1 2.3 1.9 2.0 4.1 
Co ps) 2.3 22 1.9 2.0 1.6 23 25 2.6 23 ae 3.6 
Cet 1.6 hdl el 0.9 0.9 0.7 i 1.6 ier 1.6 135 12 
Gas wetness 0.45 0.46 0.49 0.51 0.57 0.66 0.40 0.39 0.40 0.38 0.43 0.69 
Thompson Indices 
A 0.35 0.32 0.30 0.30 0.28 0.26 0.30 0.32 0.31 0.32 0.32 0.31 
B 0.87 0.68 0.67 0.66 0.65 0.63 0.61 0.74 0.72 0.74 0.72 0.54 
4 0.44 0.25 0.33 0.55 0.37 0.00 0.22 0.34 0.29 0.27 0.25 0.23 
C 1.08 1.15 1.19 1.24 1.19 1.28 1.14 1.11 1.10 1.08 1.09 129 
I 1.95 2.13 2.13 2.11 2.08 2.07 2.19 222 “22 ROD) 2.20 2.10 
F 0.58 0.63 0.63 0.62 0.59 0.60 0.65 0.68 0.68 0.68 0.68 0.66 
H 16.5 16.4 16.4 15.3 15:3 14.8 17.4 18.8 18.8 19,2 19.0 15.3 
U 1.16 1.14 1.11 1.08 1.08 1.04 Lae 1.20 1.20 122 1.20 1.06 
R Le72 1.63 1.61 1.51 [352 1.47 1.67 1.85 1.85 1.92 1.9] 1.49 
S 20.4 18.1 175 16.4 14.7 12.9 20.1 222 213 22.4 21.7 18.1 
b. Oils and condensates 

Oil Condensate 

S. Flash Stage 1 Stage 2 Stage 3 Stage 4 Stage 5 Stage 6 Stage 1 Stage 2 Stage 3 Stage 4 Stage 5 Stage 6 
Composition (wt% ) 
C,-C; 1.0 1.0 0.7 0.9 1.0 1.3 0.2 13 0.7 13 1.0 0.1 8.6 
Ceo ios il 1.0 eal 1.1 1.1 0.3 LS 1.4 2.0 1.8 0.5 9.4 
Cr 2.9 3.0 2.8 2.8 2.8 257 Ld 4.6 4.7 6.8 6.2 2.8 21.8 
Ce 4.9 pal 4.6 4.5 4.3 4.1 2.4 8.9 9.9 Foe 14.7 8.4 27.7 
CoS ee 3.8 3.5 3.4 3.1 3.0 2.3 deo 8.9 14.7 16.0 10.9 14.4 
Ciot 86.4 86.1 87.4 87.4 87.8 87.8 93.7 76.2 74.4 60.1 60.3 77.3, 18.2 
Thompson Indices 
A 0.58 0.58 0.58 0.59 0.58 0.58 0.80 0.50 0.50 0.49 0.48 0.50 0.52 
B 1.70 Lay 172 eZ 1.72 1.76 2.31 1.54 [52 1.49 1.50 1.54 1.41 
4 12) 1.23 122 1.28 125 127 1.57 1.12 1.07 Tad 1.14 1.16 0.97 
C 0.61 0.61 0.59 0.60 0.60 0.60 044 064 40.64 £40.66 0.66 0.60 0.70 
I 1.35 1.38 135 1.35 1.33 1.32 LS 1.64 1.49 1.54 1.72 1.49 1.45 
F 0.51 0.51 0.49 0.49 0.49 0.48 0.39 0.56 0.57 0.59 0.60 0.59 0.57 
H 19.8 19.7 19.3 19.1 18.8 18.4 17.2 21.8 22.8 23 23.9 24.6 20.0 
U Lea a 72 1.68 1.66 1.65 1.96 La 1.80 1.78 1.80 2.06 1.56 
R 3.21 3.15 3.17 3.11 3.05 3.06 3.49 2.81 3.50 3.54 3.06 4.04 2.84 
S 48.2 46.0 50.3 43.8 43.6 35.3 50.7 54.5 52:1 55.8 58.7 103s 41.1 


* Gas wetness = sum(C,—Cs;)/Sum(C,—Cs). 
> Thompson indices: C, (n-hexane + n heptane)/(cyclohexane + methylcyclohexane); U, cyclohexane/methylcyclopentane; A, benzene/ 
n-hexane; I, (2-+3-methylhexane)/Sum(dimethylcyclopentanes); R,n-heptane/2-methylhexane; B, toluene/n-heptane; R, n-heptane/ 
methylcyclohexane; S, n-hexane/2,2-dimethylbutane; XY, m/p xylene/n-octane; H, 100 * n-heptane/sum(heptanes) (chromatographically 


defined). 


© **Pseudo-components”’ combining all peaks in a given elution range, e.g. C5 includes all material eluting after nC, up to and 


including nCS. 
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Fig. 4. Bulk composition of oils, condensates and gases. 


of gas with a high methane content during each stage 
and the subsequent removal of the gas phases. Also 
here, the largest changes occur at the lower pressures. 

Gas wetness [o(C2—Cs)/o(C,—Cs); Table 2] is higher in 
associated gases than in free gases, as can be expected 
from a gas-liquid equilibrium. The wetness difference 
between the free and the associated gases is in agree- 
ment with the phase fractionation effect as described by 
Clayton (1991) in his gas classification plot of methane 
content vs. 6!°C of methane. 


3.2. Gas isotope ratios 


The 5!3C curves for methane through propane (Table 
3, Fig. 5) show little change in the carbon isotope ratios 
from Stage | to 5 with no significant differences between 
the free and associated gases. The free gas from Stage 6 
(no associated gas was released) is isotopically much 
heavier (12%o for methane, 3%o for ethane, 1.6%o for 
propane) compared to the gases from the previous 
stages. The heavy isotopic signature of the Stage 6 free 
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gas could be the cumulative result of a small isotopic 
fractionation (with lighter isotopes preferentially in the 
gas phase) during each step, with the remaining dis- 
solved gas becoming gradually heavier. The overall 
effect of this would be revealed during Stage 6 when the 
remaining amount of associated gas is converted into 
free gas (note that no associated gas is present at this 
stage). However, one would have expected to see a 
(small) difference between the isotope ratios in the free 
and the associated gases from the earlier stages and this 
is not the case. Unfortunately, no isotopic ratios could 
be obtained for the associated gas from Stage 5. These 
would be expected to have similarly heavy values as are 
observed in the free gas from Stage 6, as the associated 
gas in Stage 5 is the direct precursor of the free gas in 
Stage 6. 

The 85!°C values in butanes (iso and normal) and in 
CO>, and the 5!8O values in CO, show similar varia- 
tions as observed for the 6!°C values of the C,-C; 
hydrocarbons, with heaviest values in the Stage 6 gas 
and no systematic difference between the free and asso- 
ciated gases (Table 3). The 6D values for methane show 
significant scatter but without a systematic pattern. 
Unlike the other isotope ratios, 0D values are not hea- 
vier in the free gas from Stage 6. 

The data clearly indicate that phase fractionation 
does not effect the gas isotope ratios, with the possible 
exception of fractionation at a very low pressure (1 bar, 
corresponding to a depth of 10 m), which removes the 
last remaining light components from the oil — a sce- 
nario that is not very likely to occur in nature. These 
results are in contrast to the work by Carpentier et al. 
(1996), who measured a phase fractionation effect in Co— 
C4 n-alkanes and modelled that this also should be 
present in smaller molecules, even though this was not 
measured in Co—Cg n-alkanes. 


Table 3 


3.3. Light hydrocarbon fraction (Cs—Cs) 


A significant body of work on phase fractionation 
effects has been published by Thompson (1987, 1988, 
1991). Thompson’s efforts were concentrated on chan- 
ges in the light hydrocarbon fractions, which were 
described in terms of “‘aromaticity’’, ““paraffinicity” and 
“normality” via a suite of compound ratios, often 
referred to as Thompson Indices. As these concepts and 
ratios are commonly used in the geochemical literature, 
they will also be used here to describe the results of our 
differential evaporation experiment (Table 2). 

The relative amount of aromatic compounds in the 
light hydrocarbon fraction compared to _ paraffins 
(“‘aromaticity’’) is described by Thompson indices A, B 
and X (see Table 2 for explanation). The variations in 
these parameters during the differential evaporation 
experiment are generally similar (Index B — toluene/ 
nC; — is shown in Fig. 6): (1) fairly stable values in each 
fraction during Stages 1—5, (11) higher ratios in liquid 
phases than in corresponding gas phases, both at high 
pressure and under standard conditions, and (111) lower 
ratios in condensates than in corresponding oils. Only 
during the last stage, when the pressure is reduced from 
100 to 1 bar (all at reservoir temperature 165°C), do the 
ratios change significantly and the difference between oil 
and condensate increases. Other light hydrocarbon 
parameters also display an anomalous behaviour during 
the last stage (see below). As discussed above, this is 
caused by the large mass transfer in this stage as a result 
of the — geologically unrealistic — combination of high 
temperature and low pressure. 

The relative amount of paraffinic hydrocarbons (nor- 
mal and branched) compared to naphthenic compounds 
(““paraffinicity’’), as described by Thompson indices F, 
C, and I, also varies rather systematically during our 


Differential evaporation experiment: stable isotope ratios in free and associated gases* 


Associated gas 


Free gas 


S.Flash Stage 1 Stage 2 Stage3 Stage 4 Stage 5 Stage 1 Stage 2 Stage 3 Stage 4 Stage 5 Stage 6 


5/3C (per mil vs PDB) 


Methane 44.4 44.2 44.4 44.1 
Ethane 31.0 31.1 31.3 31.0 
Propane 28.8 28.6 25.7 28.8 
Isobutane 28.7 28.5 29.1 28.2 
n-Butane 28.6 29.2 29.0 28.8 
Carbon dioxide 14.2 14.5 14.8 15.0 
dD (per mil vs SMOW) 

Methane 208 230 208 222 
5/8O (per mil vs SMOW) 

Carbon dioxide 14.4 14.6 14.6 14.2 


44.4 n.m.> 44.4 44.2 44.3 44.3 44.4 32.1 
30.9 nm. = be 30.9 31.0 3 30.8 27.6 
28.9 nm. 28.8 28.9 28.5 28.6 28.8 21 
30.4 nm. 28.4 29.1 28.8 28.0 219 Pia ee) 
28.9 n.m. 28.9 29.9 29:2 28.1 28.9 28.0 
14.4 nm. 14.6 14.3 14.1 14.0 14.8 13.0 
229, nm. 229 228 221 211 221 218 
12.9 nm. 14.3 15.0 14.0 14.0 14.7 a 


* The table also contains isotope ratios measured on a single flash reservoir gas. 
n.m.; not measurable due to small sample amount. 


b 
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experiment (Index F- nC,/methylcyclohexane — is 
shown in Fig. 7): (4) fairly stable values in all fractions 
during Stages I—S, (11) higher values in gas phases than 
in corresponding liquid phases, both at high pressure 
and under standard conditions, and (111) higher values in 
condensates than in corresponding oils. Also here, the 
largest changes in the ratios occur during Stage 6. 


Thompson Index H (nC7 as% of chromatographically 
defined Cy-range) shows a more complicated pattern 
(Fig. 8): at elevated pressures and temperatures the gas 
phase has higher H values, while under standard condi- 
tions (15°C, 1 bar) liquids have higher values than cor- 
responding gas samples. Condensates have higher values 
than oils. This behaviour illustrates clearly that both the 
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character and magnitude of phase fractionation effects 
will depend on the actual temperature and pressure of 
the system (in addition to the composition and relative 
amounts of both phases). One will have to be careful in 
constructing rules that have “universal” applicability. 


Stage: 6 5 4 


The relative amount of normal compared to branched 
compounds (“normality”) is described by Thompson 
indices U, R and S (Index R-nC;/2-methyl-hexane — is 
shown in Fig. 9). As expected, gas samples are enriched 
in branched compounds (with lower boiling points), 
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Fig. 8. Paraffinicity: Thompson index H (n-heptane as % of total C; components). 
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while the indices tend to be slightly lower in oils than in 
condensates, excluding the Stage 6 samples, where the 
condensate has lower values than the oil. This varying 
behaviour again illustrates the difficulty of establishing 
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general rules for phase fractionation effects on indivi- 
dual geochemical parameters. 
Thompson (1987) introduced a plot of Index F against 
Index B to recognise evaporative fractionation (and other 
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Differential evaporation experiment: gross composition and carbon isotope ratios of whole fluids and fractions 


Oils 


Condensates 


S.Flash Stage 1 Stage 2 Stage 3 Stage 4 Stage 5 Stage 6 Stage 1 Stage 2 Stage 3 Stage 4 Stage 5 Stage 6 


Gross composition (wt.% ) 


Saturates 76 78 77 78 78 
Aromatics 20 17 20 18 17 
Polars 3 4 2 2 3 
Asphaltenes 0.4 l 12 0.9 3 
5/3C (per mil vs PDB) 
Whole fluid 27.8 28.1 28.7 28.4 28.2 
Saturates 28.8 28.6 28.6 28.9 28.8 
Aromatics 27.3 214 292° 273 21 
Polars 27.6 213 217 “215° 270 
Asphaltenes 28.0 28.1 28.2 28.5 28.2 

* Unreliable due to small sample amount. 

b n.d., not detectable. 
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80 75 85 83 80 78 76 fi) 
16 20 13 16 iy 16 21 18 

3 ) 2 2 6 =) 7 

Id 2.4 0.1 0.2 0.4 0 0 0 
28.3 28.1 28.5 284 284 27.9 n.d.> 28.2 
20.9 26.9 25:8 23.9 289° 28.0" -26.7 28.4 
2h2> 2h: 20.9" 209, 202. 25s - 268° 279 
Zbl 230) “2755: 220 279 264 274 Ziad 
28.3 28.1 28.0 280 28.4 n.d. nd. 291° 
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Fig. 11. 5'°C of whole fluids and saturated hydrocarbon fractions. 


processes). Plotting our data in a similar way (Fig. 10) 
shows relatively small differences between oils and con- 
densates, particularly after exclusion of the oil from the 
geologically unrealistic Stage 6. These results are similar 
to those reported by Thompson (1987) and Larter and 
Mills (1991) for single-step phase fractionations (see 
arrows in Fig. 10). Much larger fractionations were 
obtained by Thompson (1987) by repeatedly equilibrat- 
ing an oil with large amounts of methane with sub- 
sequent removal of the gas phase and similar large 


variations were also observed in several natural systems 
(Dzou and Hughes, 1993; Meulbroek et al., 1998). 
While a scenario of repeated addition of gas to a reser- 
voired oil may be realistic for the Gulf of Mexico, we 
believe that our relatively simple phase fractionation 
experiment without addition of excess dry gas is more 
applicable to the situation on the Norwegian con- 
tinental shelf. 

Overall, our observations on the Cs6—Cy hydrocarbons 
are qualitatively similar to those reported by others (e.g. 
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Fig. 12. Gas chromatograms of saturated hydrocarbon fractions from oils and condensates. 


Thompson, 1987; Larter and Mills, 1991) with con- 
densates having higher “paraffinicity’ and lower 
“aromaticity” and “normality” than corresponding oils. 
Quantitatively, our oil-condensate differences are smal- 
ler than those found in the other studies, primarily 
reflecting the experimental design (“‘simple’’ phase equi- 
librium without addition of dry gas) and the composi- 
tion of the starting fluid. 


3.4. Gross composition of oils and condensates 


The gross composition of the oil samples does not 
change much over the entire experiment (Table 4). In 
fact, the observed differences are generally in the same 
order as the reported uncertainty in Iatroscan data 
(approx.+1—2%; Bharati et al., 1997). Tentatively, con- 
densate samples are slightly richer in saturated hydro- 
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Table 5 
Differential evaporation experiment: GC and GCMS parameters from saturated and aromatic hydrocarbon fractions* 


Oils Condensates 


S.Flash Stage 1 Stage 2 Stage 3 Stage 4 Stage 5 Stage 6 Stage 1 Stage 2 Stage 3 Stage 4 Stage 5 Stage 6 


Pristane/nC7 0.62 0.65 0.63 0.63 063 £0.63 0.62 064 0.71 0.74 0.71 0.71 0.69 
Phytane/nCj¢ 0.43 044 044 043 O43 O47 043 O43 044 049 046 £0.47 ~~ # 0.44 
Pristane/phytane 1.54 1.49 1.58 157 1.58 1.42 1.54 1.71 i 195° 2.05 1.94 2.13 
NC 17/(NCy7 + nNC7) 0.63 0.60 062 0.62 062 0.61 0.62 0.75 0.77 0.91 0.96 0.95 0.93 
MSF (molar slope factor) 1.17 LVF 1.18 LV 1218 1.16 LJ 123 1.25 L.37 1.46 1.38 1.41 
MPII1 0.71 0.69 0.65 0.71 0.73 0.71 0.65 0.68 0.68 0.67 #£4x0.62 065 0.64 
20S 0.61 0.58 0.61 0.60 057 £0.56 O57 £40.56 0.62 0.58 0.46 0.47 n.d. 
BB 0.63 0.62 0.63 0.63 060 0.61 0.61 0.60 062 062 058 O56 ~~ nd. 
Da Ts 456 419 435 389 464 400 4.31 486 454 3.27 ~ nd. n.d. n.d. 
4R/3R 0.59 0.53 O55 0.56 0.61 0.61 0.59 047 044 0.39 0.36 £0.83 nd. 
27/29st 0.84 085 0.83 083 0.87 40.85 0.81 0.92 0.88 0.90 1.05 1.20 = nd. 
3R/27dia 0.29 0.35 0.30 0.31 0.29 0.31 0.31 0.39 0.36 0.53 0.81 1.06 nd. 
Dia/reg 1.90 E93. «2.04 2.00 221 2.26 2.200 2.38 2.49 2.29 2.70 3.20. nd. 


* Molar Slope Factor (MSF)=molar concentration ratio of two successive n-alkanes, averaged for nC,5 to nC5 (after Thompson, 
1999a,b); MPI1 =1.5*(3MP+2MP)/(P+9MP+ 1MP); P=phenanthrene, MP=methylphenanthrene; 20S=20S/20(R+ S) 140,17a- 
C>9-steranes; BB=14B,17B/(14B17B + 140,17) Cro-steranes 20(R+S); 4R/3.R=Co,-tetracyclic terpane/C3-tricyclic terpane; 27/ 
29st = Co7/Cro 50,148,17B-steranes 20(R +S); 3R/27dia = C>3-tricyclic terpane/13B,17a-C>7-diasteranes 20(R+ 8S); Dia/reg = 13B,17a- 
C,7-diasteranes 20(R + S)/14B,17B-C>7-steranes 20(R+ S); measured in m/z 259 and m/z 218, respectively. 

> n.d., not detectable. ) 
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Fig. 13. Pristane/nC,7 and pristane/phytane ratios in oils and condensates. 
carbons (and correspondingly lower in aromatic hydro- In general, our samples show less variation in gross 
carbons) during the earlier stages and gradually become composition data than those reported from other 
more like the oils. This effect 1s in accordance with a experiments; e.g. Larter and Mills (1991) found up to 
higher “‘dissolving’” power of the gas phase for long- 10% enrichment of saturated hydrocarbons in their 


chained (>C,,) paraffins at high pressures (see below). condensate samples. 
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3.5. Carbon isotope ratios of whole oils/condensates and 
fractions 


5'3C values for whole oils/condensates and _ their 
fractions (Table 4) show relatively small changes 
throughout the differential evaporation experiment. As 
a rule, the corresponding 5!°C values for oils and con- 
densates are within 0.5%o of each other, which is only 
slightly more than the analytical uncertainty, particu- 
larly with the small sample size that is the case for sev- 
eral of the condensate fractions. 

The 6!°C values for the saturated hydrocarbons, 
which are probably the most reliable due to the larger 
size of this fraction, show an interesting pattern with 
lighter ratios in the condensates compared to the oils in 
Stages | and 2 and heavier condensates during the latter 
stages (Fig. 11). This pattern is in agreement with the 
change in the composition of the saturated hydrocarbon 
fraction (see below) with relatively more (isotopically 
light) long-chain alkanes in the condensate samples 
from Stages | and 2 and less in the later condensate 
samples. Note that the observed variations are small 
and in some cases come close to the uncertainty in the 
analytical data (normally <0.1%o, but can be larger in 
small samples). 

Differences of max 0.5%o between the isotopic ratios 
of whole fluids and fractions of oils and condensates, 
as observed in our experiments, would be considered 
as not significant in most correlation studies. However, 
Dzou and Hughes (1993) found isotopic shifts of up to 
1.5%o due to phase fractionation. These were caused 
by the very waxy nature of their oils and the resulting 
loss of isotopically light long-chain n-alkanes in the 
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condensates, but if not recognised as phase effects, 
these differences could easily have been interpreted as 
indicating different sources for the oils and con- 
densates. 


3.6. Cy5+ saturated hydrocarbons 


The gas chromatograms of the saturated fractions 
from all oil samples (Fig. 12) are very similar and show 
a relatively large contribution of long-chain n-alkanes 
(>C,>5). In the condensates, the amount of long-chain n- 
alkanes is slightly lower than in the oils during the first 
two stages at high pressures and this decreases rapidly at 
pressures below 300 bar (Stage 3 and onward). This 
clearly reflects a limited “‘solvent power” of the gas 
phase for longer chain molecules, as well as a decrease 
of this solvent power at lower pressures. 

The lower solvent power of the gas phase for larger 
molecules at low pressures is also recognisable in the 
various parameters that are derived from the gas chro- 
matograms (Table 5). Pristane/phytane ratios are 
slightly higher in condensates than in oils at high pres- 
sures, but this difference increases at lower pressures 
(Fig. 13). Note, that the Pr/Ph ratio of the oil does not 
decrease noticeably during the experiment, due to the 
small absolute quantities of pristane and phytane that 
are evaporated from the oil. As can be expected, the 
C17/Ci7+ C7 ratio shows higher values in condensates 
compared to oils, with an increase at lower pressures 
(Fig. 14). The stability of this ratio in the oil samples 
again illustrates the small amounts of these (relatively 
heavy) compounds that are evaporated. 
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Fig. 14. Molar Slope Factor (MSF) and C,7/(C,7+ C27) ratio in oils and condensates. 
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Phytane/nC 1g ratios are similar in oils and condensates 
(Table 5) and do not vary during the experiment. This is 
in line with the similar retention times (strongly depen- 
dent on boiling point) for both compounds. Somewhat 
surprisingly in view of the similar relative retention 
behaviour as phytane and nCjg, pristane/nC;7 ratios are 
slightly but consistently higher in condensates than in oils 
for Stages 2 to 6 (Fig. 13). This illustrates the uncertainties 


23-3R 


| 24-4R 


Javenennne 
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in extrapolating boiling point differences (retention times) 
from standard conditions (1 bar, 15°C) to significantly 
higher temperatures and pressures. 

Changes in the overall appearance of gas chromato- 
grams as described above have also been reported by 
Thompson (1987, 1988), Larter and Mills (1991) and 
Knudsen and Meisingset (1991). Larter and Mills have 
also evaluated the Pr/Ph, Pr/nC,7 and Ph/nCjg ratios, 


m/z 191 
terpanes 


m/z 217 
steranes 


Fig. 15. Biomarker distributions (m/z 191 and m/z 217) in saturated hydrocarbon fractions from the starting reservoir oil. 23- 
3R =C,3-tricyclic terpane; 24-4R = C,,-tetracyclic terpane; T, = 180-22,29,30-trisnorneohopane; 30 dia = 15a0(Me)-17a-27-norhopane; 


3008 = 170,21B8-hopane; 27 BD = 13B,17a-diacholestane; 29BD 


= 24-ethyl-13B,17a-dicholestane; 29 BB = 24-ethyl-148,17B-cholestane. 
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but unlike our results, they did not observe clear differ- 
ences between oils and condensates. Higher Pr/Ph 
values for condensates compared to oils were reported 
by Dzou and Hughes (1993). In a study of samples from 
Vermilion 39, Curiale and Bromley (1996) identify two 
groups of samples and attribute the differences between 
the groups to migration fractionation processes. These 
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differences include variations in the total gas chromato- 
grams, similar to those found in our experiments, and 
significantly higher Pr/Ph ratios (shift from 2—2.5 to 3— 
4.5) in fractionated samples. 

In work on oils and condensates from the Gulf of 
Mexico, Thompson (2000a,b) introduced the so-called 
Molar Slope Factor (MSF) as a measure of the decrease 
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Fig. 16. 148,178/(14B,17B + 140,170) Coo steranes (bb) and C57/Cr9 148,17B-steranes (27/29) in oils and condensates. 


3 2 4 Reservoir 


ee ee ee fli 


Stage: 6 5 4 
3.5 
3 
225 
s 
G 2 
tee 
0 
8 1.5 
mo) 
oc | 
~* 
0.5 
0 
0 100 200 


300 400 500 


Pressure (bar) 


* 3R/dia cond * 3R/dia oil -* dia/bb cond * dia/bb oil 


Fig. 17. C3-tricyclic terpane/C>7-diasteranes (3R/dia) and C>7-diasteranes/C>7-148,17B-steranes (dia/bb) in oils and condensates. 


1436 G.W. van Graas et al. | Organic Geochemistry 31 (2000) 1419-1439 


in n-alkane concentrations with increasing chain lengths 
that normally occurs in mature marine oils.4 According 
to Thompson’s results “normal” oils have slope factors 
lower than 1.2, with thermal condensates having slope 
factors of approximately 1.2 and evaporative con- 
densates having substantially higher values. Our data 
are in agreement with this (Fig. 14) with constant slope 


OILS 


495 bar 


=e 


390 bar 


factors in the oils (1.16—1.17) and significantly higher 
values in the corresponding condensates (1.23—1.46), 
particularly at lower pressures (< 300 bar). This sup- 
ports Thompson’s suggestion that the MSF can be a 
useful tool in the classification of oils and condensates. 

The reservoir fluid used in our experiment has a rela- 
tively high maturity (deep end of the oil generation 
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Fig. 18. Gas chromatograms of aromatic hydrocarbon fractions from oils and condensates. 


+ According to Kissin (1987) the logarithmic molar concentrations of n-alkanes in petroleum are linearly related to their carbon 
number. In such a plot, the slope factor would express the angle of the straight line. As used in this paper, the MSF is equal to the 
molar concentration ratio of any pair of adjacent alkanes: the MSF=[n]/[n+ 1]. By averaging individual values, slope factors can be 


calculated for different carbon number ranges (see Table 5). 
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Fig. 19. Methylphenanthrene index (MPI1) in oils and condensates. 


window), which is clearly reflected in the biomarker 
distributions (high diasteranes, 7,/7,, and diahopane; 
see Fig. 15). The quality of the GCMS data for the 
saturated hydrocarbon fractions is generally good with 
the exception of the condensate samples from Stages 4-6 
which have very low biomarker concentrations and 
correspondingly poor signal/noise ratios, particularly in 
case of the triterpanes. Only a few biomarker com- 
pounds could be detected in those samples (none in the 
condensate from Stage 6). 

In general, smaller molecules are enriched in con- 
densates compared to oils and individual biomarker 
parameters are affected accordingly. Ratios calculated 
from molecules of different size or from isomers with a 
different retention behaviour, such as C33 tricyclic ter- 
pane/C>,7 diasteranes, Cs7/Co9 BB-steranes or dia/BB Co7 
steranes, show fairly similar values in oils and condensates 
at high pressures but increasingly higher values in con- 
densates at lower pressures while the values in oils 
remain unchanged (Figs. 16 and 17). In contrast, a ratio 
composed of more similar molecules, like BB/(a« + BB) 
Cro steranes (Fig. 16), has similar values in oils and 
condensates throughout the entire experiment. Depend- 
ing on the pressure, several biomarker ratios will show 
differences between coeval oils and condensates that 
could easily lead to erroneous conclusions in correlation 
studies if not recognised as being induced by phase 
fractionation. This conclusion was also reached by 


Curiale and Bromley (1996), who observed variations in 
biomarker distributions very similar to those seen in our 
samples. 


3.7. Cy5+ aromatic hydrocarbons 


The gas chromatograms of the aromatic hydrocarbon 
fractions from all oil samples (Fig. 18) are rather simi- 
lar, as is the case for the saturated hydrocarbon frac- 
tions. The differences that can be observed, are probably 
mainly due to variations introduced during the solvent 
evaporation steps where a variable amount of the low- 
boiling compounds (C,- and C>-naphthalenes) can be 
lost. More variation is seen amongst the condensate 
samples. In general, condensates differ from oils in hav- 
ing a lower Unresolved Complex Mixture (UCM) hump 
at high retention times. As in the case of the saturated 
hydrocarbons, this difference is more pronounced at 
lower pressures, when the solvent power of the gas 
phase for these relatively large aromatic compounds 1s 
greatly reduced. 

Inspection of the behaviour of an important aromatic 
maturity parameter, the methyl- phenanthrene index 
(MPI1; Radke, 1988), shows generally slightly higher 
values in oils than in condensates (Fig. 19). However, 
the differences are small and of the same magnitude as 
the analytical uncertainty. This is due to the way in 
which the MPI1 is calculated with the compounds in the 
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numerator (3- and 2-methylphenanthrenes) having an 
intermediate GC-retention time behaviour compared to 
the compounds in the denominator (phenanthrene and 
9- and l-methylphenanthrenes). Insignificant variations 
in MPI1 values were also reported by Larter and Mills 
(1991) and Dzou and Hughes (1993). In contrast, 
Thompson (1991) did report significant differences 
between oils and condensates as a result of phase frac- 
tionation, but in his case the fractionation conditions 
were much more severe. 

Due to the relatively high maturity of the reservoir 
fluid that was used in our experiments none of the sam- 
ples contains detectable amounts of C5¢.-Cr9 mono- or 
triaromatic steroid hydrocarbons. 


4. Conclusions 


In a six-stage differential evaporation experiment only 
small variations were observed in the composition of the 
oil samples. This is largely due to the composition of the 
starting fluid (gas—oil ratio 400 Sm*/Sm2?), giving a 
quantitative dominance of the liquid phase throughout 
the experiment. 

The compositions of the evolved free and associated 
gases vary as can be expected, with the associated gas 
being wetter than the free gas and becoming increasingly 
wet at lower pressures. Carbon isotope ratios of C,;—Cy, 
hydrocarbons show no significant differences between 
the free and associated gases. 

Changes observed between the light hydrocarbon 
(Ce—Co) fractions of the oils and condensates are similar 
to those reported by other workers (e.g. Thompson, 
1987) with oils having lower “normality” and higher 
“aromaticity” and “‘paraffinicity” than the correspond- 
ing condensates. Careful evaluation of the data shows 
that not only the magnitude of phase fractionation 
effects is dependent on temperature and pressure, but 
also that the direction of fractionation effects (i.e. if 
parameters increase or decrease) may change when 
going from high temperature and pressure to standard 
conditions (15°C, 1 bar). 

Evaluation of the saturated and aromatic hydro- 
carbon fractions in condensates shows that the gas 
phase has a good solvent power for large (> C19) mole- 
cules at high pressures (>300 bar), which is con- 
siderably reduced at lower pressures. Differences 
between oils and condensates for individual parameters 
are very much as can be expected from molecular size 
(boiling point, GC-retention time) of the compounds 
involved, with smaller molecules being preferentially 
enriched in the condensates. Phase fractionation effects 
on several geochemical parameters can be sufficiently 
large to interfere with a “simple” correlation between 
oils, condensates and source rocks. 


In general, it will be very difficult to set up detailed 
guidelines for how individual geochemical parameters 
will change due to phase fractionation, since these will be 
controlled by a complex set of parameters: (1) pressure 
and temperature of fractionation, (11) composition of the 
fluids involved, e.g. API gravity of the oil, composition of 
the gas, and (11) relative amounts of oil and gas in the 
system (gas—oil ratio of a single-phase petroleum, 
amount of gas added to an oil-filled reservoir). Different 
petroleum fractions (and their derived parameters) will 
be affected most strongly depending on the specific 
conditions under which phase fractionation occurs. 

In summary, the results obtained from our experiments 
further support the importance of phase fractionation as 
an important process to consider when interpreting 
geochemical data for both condensates and oils. 
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